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Abstract 

We consider the classical XY model (0(2) nonlinear a- model) on a class of lattices 
with the (fractal) dimensions 1 < D < 2. The Berezinskii's harmonic approximation 
suggests that the model undergoes a phase transition in which the low temperature 
phase is characterized by stretched exponential decay of correlations. We prove an 
exponentially decaying upper bound for the two-point correlation functions at non- 
zero temperatures, thus excluding the possibility of such a phase transition. 

The classical XY model (or, equivalently, 0(2) nonlinear a-model) has been a subject of con- 
siderable interest in the contexts of statistical physics and relativistic field theory. The model 
exhibits a standard ferromagnetic phase transition in dimensions d > 2, while it undergoes 
an exotic phase transition called the Berezinskii-Kosterlitz-Thouless transition [1, 2, 3] in 
d = 2. The existence of a phase transition in two dimensions is in a remarkable contrast 
with the 0(n)-rotator models (a-models) with n > 3 such as the classical Heisenberg model, 
which are conventionally believed to be asymptotically free and have no phase transitions in 
two dimensions [4] . In one dimension, general arguments guarantee that the XY model (like 
all the other short range spin systems) has no phase transitions. 

It is now common to regard the dimension dasa continuous parameter, having in mind 
lattices with fractal structures for example. Then a natural question is whether the classical 
XY model exhibits a (finite temperature) phase transition in the intermediate dimensions 
1 < d < 2. Although there have been no publications directly devoted to this problem as 
far as we know, we believe that the problem is important and is worth settling. 

To illustrate that the problem is nontrivial, denote by T c (d) the critical temperature (at 
which the susceptibility diverges) of the classical XY model on the ci-dimension hypercubic 
lattice. The rigorously established facts are that T c (d) is finite for d = 2,3,4, ■■ •, and is 
vanishing for d — 1. If one naively regards T c (d) as a function of the continuous parameter 
d, then the most natural "interpolation" may be that T c (d) continuously decreases below 
d = 2, and vanishes continuously at d = 1 (or some dimension 1 < d < 2). A much stronger 
support for this naive guess comes from the observation that a straightforward extension of 
the Berezinskii's harmonic approximation [1] (which correctly predicts the existence of the 
Kosterlitz-Thouless transition in d = 2) indicates that, in 1 < d < 2, there can be an exotic 
low temperature phase characterized by stretched exponential decay of correlations. One 
who came across with this argument might conjecture the existence of a phase transition in 
1 < d < 2. 

1 Published in Phys. Rev. Lett. 74, 3916 (1995). 
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In the present Letter, we address the above problem and provide a conclusive result. 
We prove that, in the classical XY model defined on a lattice with the (fractal) dimension 
1 < D < 2, the two-point correlation function decays exponentially at any finite temperature. 
Our (fractal) dimension D coincides with the euclidean dimensions for regular lattices. This 
result excludes the possibility of the new type of phase transition mentioned above [5]. It 
also leads us to a somehow unexpected conclusion that T C (D), as a function of continuous 
dimension D, drops discontinuously [6] to zero at D — 2. 

The proof is carried out in two steps. First we use the McBryan-Spencer method [7] 
to prove an upper bound for the two-point correlation function in terms of a stretched 
exponentially decaying function. This result automatically extends to a very general class 
of models with continuous symmetries, and strengthens the well-known result that these 
systems have no long range order in d < 2 [8, 9, 10]. Next we apply the Simon's inequality 
and his argument on the decay of correlation [11] to conclude that the correlation decays 
exponentially in the classical XY model. 

Definitions: We consider a general connected lattice A = (A s ,A b ), where A s is a set of 
countably infinite sites x,y, . . ., and A b is a set of bonds, i.e., pairs of sites (x, y), (u, v ), . . . 
through which spin variables directly interact. We define the classical XY model on A. The 
Hamiltonian (action) is given by 

H = - s «- s « = - H cos(9 u -9 v ), (1) 

(u,v)eA b (u,v)eA b 

where s x = (cos 9 X , sin 9 X ) is the spin variable at site x G A s . The two-point correlation 
function at inverse temperature (3 > is 

(s x ■ Sy) p := Zf x (n de)j e~ m cos(9 x - 6 y ) (2) 

with the partition function 

^:=(ll£^)e-^. (3) 

The expectation value (2) is to be interpreted as the thermodynamic limit of the correspond- 
ing finite volume quantities. 

Berezinskii's argument: Before describing our results, we shall briefly review the Berezin- 
skii's argument and discuss its extension. The essence of the Berezinskii's harmonic approx- 
imation [1, 7, 12] is to make the model into a Gaussian model by replacing cos(6> u — 6 V ) in 
the Hamiltonian (1) by [1 — (9 U — 9 v ) 2 /2] and replacing the range of integrals in (2) and (3) 
to (— oo, oo). A naive expectation is that the approximation gives sensible results for large (3 



2 



where local difference of fields (9 U — 9 V ) is small [13]. After making the above replacements 
(and suitably eliminating contributions of a zero- mode), the correlation function is evaluated 

as 



la e \ Gauss 



(rLHoO e- pn cos(9 x - e y ) 



= exp[-((9 x - 9 y ) 2 )f^/2] = explp-'C^y)^}, (4) 

where TL = J2( u ,v)eA b (9u — 9 v ) 2 /2 = — J2ueA B #«(A6>) u /2 is the approximate Gaussian Hamil- 
tonian. Here the lattice Laplacian is defined as (A/) z = J2 y eA s :( y ,z)eA h (fy~fz)- The quantity 
C(x,y) = —f3((9 x — 9 y ) 2 )^ auss can be characterized as follows. For fixed x, y, let ip z be the 
solution of the lattice Poisson equation 

- (Ay?) 2 = S XtZ - 8 VtZ . (5) 

Then C(x,y) is equal to the "potential difference" C(x,y) — ip y — ip x . By examining the 
solution of (5) in a general ci-dimensional lattice, we find that C(x, y) behaves for large \x — y\ 

as 

' A(d)\x - y\ 2 ' d - B(d) if 2 < d 
C(x,y)~l - A(d) log \x - y\ if d = 2 , (6) 

k -A{d)\x-y\ 2 - d if 1 <d<2 

where A(d) and B(d) are strictly positive constants which depend only on the lattice struc- 
ture. The asymptotic behavior (6) for nonintegral d was obtained by extrapolating from 
those for integral d. On a fractal lattice, one should interpret d in (6) as a suitable "fractal 
dimension" [14]. 

The approximation (4) along with (6) implies that the asymptotic behavior as \x — y\ | oo 
of the correlation functions for large f3 (low temperatures) is given by (s x • s^)^ ~ °"(/5) 2 > 
for d > 2, (s x • s y ) a ~ \x — y^W) f or d = 2, and (s x ■ s y ) ~ ~ exp[— \x — y|/£(/3)] for d — 1. 
Here (the approximations for) the order parameter a(j3), the critical exponent r)(/3), and 
the correlation length £((3) are given by a((3) = exp[-S(d)/(4/3)], r]((3) = A(2)/(2(3), and 
£,((3) = 2/3/A(l). These estimates for correlations (including the /5-dependence of a(/3), 
7](/3), and £(/?)) recover the known (or expected) behavior in a semi-qualitative manner. It 
is remarkable that the simple approximation yields such strong results. 

In dimensions 1 < d < 2, the same approximation yields the asymptotic behavior 

(s x -s y ) p ^exp[-a((3)\x-y\ 2 - d ], (7) 

with a positive function a((3) of < (3 < oo, which is usually referred to as a stretched 
exponential decay. Since the correlation decays exponentially for sufficiently small (3, this 
observation suggests that the classical XY model undergoes a phase transition in 1 < d < 2. 
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It is interesting that, as we will show in Lemma, the McBryan-Spencer argument leads us 
rigorously to an upper bound for the correlation decaying by a similar stretched exponential 
law. 



Main results: We define the "sphere" S n (x) of radius n centered at x G A s by 

S n (x) := {y E A s |dist(x,y) = n}. 




Here dist(x, y) is the graph-theoretic distance between the sites x, y, which is defined as the 
minimum number of bonds in Ab that one needs to connect x and y. We assume that there 
exists a "(fractal) dimension" D of the lattice [15] such that the number of sites in S n (x) is 
bounded uniformly from above as 



with a positive constant C. It is obvious that the dimension D coincides with the euclidean 
dimension if A is a regular lattice. 

Then our main result is the following exponentially decaying upper bound for the two- 
point correlation. 

Theorem — If the dimension D satisfies 1 < D < 2, the two-point correlation function is 
bounded as 



for any < f3 < oo, where m(j3) > and R(/3) > are functions of f3 [16]. 

In the proof of the theorem, we make use of the following lemma. The lemma states an 
upper bound for the two-point correlation in terms of a stretched exponentially decaying 
function like the one obtained from the Berezinskii's approximation. 

We note that the following lemma can be trivially extended to cover more general class 
of systems with a global continuous symmetry. The examples include 0(n) (n > 2) rotators 
(nonlinear a- model) [7], quantum Heisenberg models [17], and the Hubbard model [18] for 
itinerant electrons. The main theorem, on the other hand, can be proved only in the classical 
XY model for the technical reason that the Simon inequality is known only for one- and two- 
component spin systems. We suspect that the statement of the theorem is valid for the 
larger class of models. 

Lemma — For 1 < D < 2, we have 




(9) 



|(s x • Sy) p \ < exp[-m(/3){dist(x,y) - R{(3)}\ 



(10) 




(11) 



for any < (3 < oo, where the function /(/3) is given by 

f{(5) = max [q - [3C 2 (cosh q- l)/(2 - D)} > 0, 
and is decreasing in (5. 



(12) 
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Proof of Lemma: We fix two sites x, y throughout the proof. Following McBryan and 
Spencer [7], we make the complex transformations 



g -> 9. + ic 



(13) 



in the numerator of (2), where <p z are real numbers to be determined. This means that we 
deform the path of integration and use the periodicity of the cosine to cancel contributions 
of the lateral contours. The transformation combined with \e lz \ = 1 for real z yields 



(s x • s^l < e-^-^Zj 1 (J] j*_je)j exp 



(3 E cos (#« - Ov) cosh(0 u - <f> v ) 
(u,v)eA b 



< e -(4>x-<t>y) 



exp 



P E ( c osh(0 u - (j) v ) - 1) 
. (u,v)eA h 



(14) 



Following the idea of Picco [19], we choose {<fi z } as 

4> z 



' q(N 2 ~ D - 1) iiz = x 

q(N 2 - D - dist(x, z) 2 ~ D ) if 1 < dist(x, z) < N 
if dist(x, z) > N 



(15) 



where N = dist(x, y), and q is a positive constant to be determined. From the choice (15), 
we have that 

0* - 2 = (16) 

for z such that dist(x, z) — 1, and that 



\4>u ~(f>v\<qn <q 



(17) 



for u, v satisfying n = dist(x, u) = dist(rr, v) — 1 ^ 0. 

The summation in the exponential in the right-hand side of (14) can be evaluated as 

E (cosh(0 u - <j> v ) - 1) 
(u,v)eA b 

N-l 

= E E E [cosh(0 u - (f) v ) - 1} 

n=l u:dist(x,u)=n v:(u,v)EA b 

dist(a;,i;)>dist(a;,«) 

<E E E (coshg-l)g- 2 (0 M -0„) 2 

n=l u:dist(x,u)=n v:(u,v)£A b 

dist(a;,u)>dist(a;,«) 
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N-1 

<C(coshg-l)£ J2 n2 ~ 2D 

n=l u:dist(x,u)=n 
N-1 

<C 2 (coshg- 1)^71^ xn 2 - 2D 

n=l 

< C 2 (coshg- l)(iV 2 ^ -l)/(2-D), (18) 

where we have used (16) and the choice (15) of {<p z } to rewrite the summation at the first 
line. The succeeding three bounds have been obtained by using (17) and (9). By substituting 
(18) and (15) into the main bound (14), we get a bound of the form (11). Finally we optimize 
the bound by choosing the constant q according to (12). 



Proof of Theorem: We assume 1 < D < 2. Recall the Simon inequality [11] for the 
classical XY model 

(19) 

uev,w£v 
(u,v)eA b 

where the finite set V is chosen so that x G V, and y ^ V. We set V = Sr(x) with some 
R < dist(x, y). By using (9) and the bound (11) of the lemma, (19) implies 

(s.-^^e-^^.s,,^ (20) 

with 

e -m(m = C 2p R D-l exp ^_j^ R 2-D _ x) ] ; (21) 

where the site v± gives the maximum value of (s^ • s^)^ in the summation in the right-hand 
side of (19). Noting that the right-hand side of (21) is less than 1 for sufficiently large R, 
we let R(j3) be the value of R which maximizes the "mass" m(/3) in the left-hand side of 
(21). We note that, for such R(f3) and for any given distance r = dist(x,y), there exist 
non-negative integers I and r' satisfying r = R(j3)£ + r' and r' < R{f3). (For r < R((3), we 
set r' — r and £ = 0.) By setting R = R(/3), we can apply the Simon inequality (19) to the 
right-hand side of (20) repeatedly at least [I — 1) more times. We thus obtain the desired 
exponentially decaying bound [20] as 

(Sx ' Sy) p < e- m ^ e (s ve ■ s y ) p < e - m ^- r '\ (22) 



Note added 2 : After the publication of the paper, we learned that the XY model on a 
class of finitely ramified fractals is studied in [21], which is one of the pioneering papers in 
"physics on fractals." It is a pleasure to thank Deepak Dhar and Tetsuya Hattori for useful 
corresp ondences . 

2 Phys. Rev. Lett. 75, 984 (1995). 
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